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Abstract 

Discrete interaction models for the classical harmonic oscillator are 
used for introducing new mathematical generalizations in the usual 
continuum formalism. Generalized discrete hyperbolic and trigono- 
metric functions are the exact solutions to generalized (varying fre- 
quency) discrete classical harmonic oscillators. Energy and angular 
momentum conservation and the preservation of the Poisson bracket 
structure with discrete interactions is discussed. 



1 Introduction 

Discrete to continuum transitions in the classical harmonic oscillator have 
been studied [1] with a simple model of a discrete harmonic oscillator where 
the continuous interaction is replaced by successive exchanges of momentum- 
energy packets, classical mimics of the interaction quanta. The subject is 
retaken here with an eye on mathematical generalizations of the usual for- 
malism as a first step towards future connections with a proper quantum 
formalism. 

Discrete to continuum transitions is a research subject in mathematics, 
physics and engineering; see, for example [2] and the references there in. 
Discretization in physics is an old and very often used tool. Generally, it is 
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a lattice structure [3] taken as an approximation of supposedly continuous 
physical systems whose continuity is then retrieved by taking a limit with 
the lattice size going arbitrarily close to zero. In such a context, discretiza- 
tion is just a mathematical artifact for overcoming computational difficulties 
with continuous systems. In contradistinction, the idea behind the naive dis- 
crete interaction models considered here is of a genuinely quantum (in the 
sense of discrete) world. Increments like At,-, Ar}, Apj, of time, position and 
momentum, respectively, although extremely small, are fixed and essentially 
non-null and can be taken as null only in an approximate sense. Nonetheless, 
the standard picture of a continuum interaction is still retrieved - otherwise 
these discrete models would be completely irrelevant- as the number of in- 
teractions grows larger, n » 1. Actually there is an interval 



for the integer n on which the continuum interaction picture is a good effec- 
tive approximation; outside it the interaction discreteness cannot be ignored. 
uj n is a system characteristic frequency and At n is the time interval between 
the n — th and the (n + 1) — th interaction. So, in its time evolution a sys- 
tem goes through three well defined phases which may present very distinct 
characteristics [1]. 

That the continuum is retrieved in the limit of zero latticing in the stan- 
dard lattice calculation is, of course, no surprising since it was tailored for 
that. Remarkable, in the case of the genuine discrete interactions, is the nat- 
ural and unavoidable emergence of an effective continuum interaction from 
fixed increments as their number n grows larger, even though there are no 
' ' forces" in between the exchange of ' ' quanta" , hence no potential energy 
and no Hamiltonian. This discrete-to-continuum transition for large values of 
n reminds the Bohr Correspondence Principle, as discussed in [1]. The limit 
of zero increment is never taken. The continuum here is the approximate 
picture. Also remarkable is the new (not contained in continuum interaction 
picture) physics outside the good-effective-approximation range for either n 
too small (ultra short distances) or too large (ultra large distances) in the 
sense of (1). From this point of view, the exact solutions of discretely inter- 
acting systems become much more interesting on the light that the standard 
continuum solutions are incomplete information [4] much in the same way 
as the quantum solutions are to the classical (the approximate) ones. Dis- 
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>> n » 1, 
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crete or continuum fundamental physical interactions is to be experimentally 
decided: Discreteness effects for n outside the interval (f) can be exper- 
imentally observed. There are, perhaps, signs that they have already [f] 
been observed. The point in [1] was to show the continuum solution of the 
standard harmonic oscillator emerging as first-order approximation in the 
large-n limit of solutions of the discrete harmonic oscillator. Discrete exact 
solutions of discrete harmonic oscillator, for all values of n, are the subject 
of this paper. It is structured in the following way. After briefly review- 
ing the idea of classical discrete interactions in Section II, we treat the time 
evolution of a discrete harmonic oscillator as a generalized discrete rotation 
in phase space introducing the generalized discrete trigonometric functions 
as its exact solutions; discrete hyperbolic harmonic oscillator (an inverted 
harmonic potential) is also considered. In Section III we deal with the sim- 
plest case of constant time interval and its transition, in the large n limit, 
to the standard continuous harmonic oscillator. Angular momentum and en- 
ergy conservation and the preservation of the Poison bracket structures are 
discussed in Section IV where the validity of Eq. (1) is established. Finally, 
with the Conclusions in Section V, we discuss the various possible alternative 
physical interpretations. 



2 Discrete interactions 

A non-relativistic regime, which is appropriate for dealing with harmonic 
oscillator, is assumed. So, we will consider the relative movement of the 
reduced mass m of a non-relativistic two body-system. At the initial time 
to, with initial position r and momentum p , it freely moves on a rectilinear 
trajectory until its, after then, first interaction at 

ti = t + At , 

— * 

n = r + Af = f H At , 

m 

that changes its momentum to 

pi = po + Ap . 
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At the n — th interaction, 

n-1 

t„ = to + E At r ( 2 ) 

3=0 

n-1 

Pn = Pn-1 + Ap n _i = p + J2 A Pj> ( 3 ) 

3=0 

p 1 U ^ 

r n = r n _i + Af n _i = f n _i + -^At n _i = f H V PjAij. (4) 

m m j^o 

The position vector, continually changing, describes a polygonal trajectory. 
There is a sudden change of momentum at each interaction point, a vertex 
of the polygon. Momentum is clearly a discrete parameter, or at least it 
changes discretely. The continuous variables time and position enter in the 
description of the motion as if they were also discrete parameters only be- 
cause the interaction events are our reference points for time counting; this is 
not a lattice calculation. The system time evolution in phase-space consist of 
changes of position produced along a non-null lapse of time but with a fixed 
momentum, and changes of momentum produced in a null lapse of time but 
at a fixed position. This is the essence of Eqs. (2-4), of discrete interactions. 
Its physical content is that any allowed apparently continuum path in phase 
space is, actually, just a first order approximation of these discrete change of 
momentum followed by changes of positions. This is a point of contact with 
Quantum Field Theories: Discreteness in the interactions through an ex- 
change of (quantum) discrete and defined amount of momentum and energy. 
One cannot talk of force or acceleration, just of sudden change of momen- 
tum or velocity. Force and acceleration, as well as all physical concepts that 
depend on them, become effective concepts, valid only in appropriate limits, 
for which the picture of a continuum interaction is a good approximation. 

2.1 Discrete Oscillator 

The Eqs. (2-4) are generic relations valid for any non-relativistic discretely 
interacting system. Let us now make Eq. (3) specific to a discrete harmonic 
oscillator DHO with the following phenomenological assumption about the 
interaction 

Ap n _i = -mw^ViAtn-i, (5) 
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which can be taken as a definition of a DHO. Eqs. (4,5) restrict the DHO 
to the plane determined by r*o and po, a s i n an y central force motion. These 
evolution equations can be put in a very convenient matrix form 



1 At n _i/m\ / f n _i 

Pn) \-17lLJ 2 At n - 1 1 / \p n -l 



(6) 



For a dimensional homogenization in the phase space we make a change of 
variable 

r' n = v^ruJ r n p' n = p n / y^muj, (7) 
so that Eq. (6), with = uAt n -i, becomes 

where we have introduced the finite discrete transformation G(9 n -\). From 
now on we will drop the primes on the new variables. Then 

G{6j) is a quasi-unimodular transformation 

det G(0j) = 1 + B) « 1, (10) 

since #j << 1, 

G(^- 1 = (l + ^)- 1 G'(-^. (11) 

The semblance of G(6j) of Eq. (8) with an infinitesimal (continuous) rotation 
suggests the introduction of generalized discrete trigonometric functions 

( cosg9 ^\ sing T\)--=ii G ^ ( i2 ) 

\— smg#(n) cosg#(n) / *-J- Q 



so that 



Therefore, 



f, n \ _ ( cosg6>(n) sing#(n) \ / r 
Pn) ~ V-smg0(n) cosg6»(n) / \p 



(13) 



n— 1 n— 1 n— 1 

cosg0(n) = 1 - £ + E ^ E •••** + 

jl<?2 il<---<J4 jl<...<j6 
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n— 1 n— 1 n— 1 n— 1 

fflng0(n) = £ E + E E • ••• 

j ii<h<jz h<—<36 h<-<37 

These generalized discrete trigonometric functions, according to Eq. (13), 
constitute the exact solutions of the discrete harmonic oscillator. They are 
both finite series and can be more compactly expressed with 

n— 1 

©?>:= E "j: (>,, ■■■<>,, 0<k<n (14) 
o<n<...<j fc 



0i o) = 1, and = if k > 



n 



so that, 



[n/2] 

co8g0(n)= E(-l) S ©i 2s) , (15) 

s=0 
[n/2] 

sing0(n):= E (~^ S+1) , (16) 
where [n/2] is the largest integer in n/2. Or yet, we can just use 

i n— 1 rt— 1 

co8g0(n) = -(11(1 +^') + 11(1 " (17) 

Z j=0 j=0 
i n— 1 n— 1 

sing0(n) = _(JI(1 + ^) - 11(1-^))- (18) 

z< j=0 j=0 
with 2 2 = — 1 and the convention that 

ni 

}[ := 1 if n2 > ni, for any Xj. 

j=m 

Then 

n-1 

cosg^(n) ± i sing#(n) = J] (1 ± i^), (19) 

3=0 

which is a generalized discrete version of de Moivre's formula. It leads to 
cosg#(n + 1) ±ismg9(n + 1) = (1 ± i9 n )[cosg9(n) ± ising9(n)], (20) 



Table 1: The generalized discrete trigonometric functions for the first five 
values of n 



cosg^(O) = 


1 




sing^(O) = 









cosg#(l) = 


1 




sing0(l) = 


00 






cosg#(2) = 


1 


— 9i9 


sing#(2) = 


00 + 01 






cosg#(3) = 


1 


— (dido + 0002 + 0102) 


sing#(3) = 


0o + 0i 


+ 02 


— 0Q0 102 


cosg#(4) = 


1 


- (e o {O 1 + e 2 + 3 )+0 1 (O 2 + 


sing#(4) = 


0o + 0i 


+ 02 


+ 03- 


+03 


H 


- 020?>) + +#0$1$2#3 


— 0q(0\02 - 


- 0103 + 


0203) 


— 010203 



f n ±ip n = [cosg6>(n) ± ising0(n)](f o ± ip ), (21) 

making explicit the change-of-phase character of the harmonic oscillator mo- 
tion in phase space. 

Table I lists the generalized discrete trigonometric functions for the first 
five values of n. 



2.2 Some properties 

This notation for the generalized discrete trigonometric functions, introduced 
in Eq. (8), explores the parallelism with the continuous formalism, like in 
the symmetry 

cosg(-0(n 2 )) = cosg6»(n 2 ) 

sing(-6>(n 2 )) = -sing6>(n 2 ), 

which corresponds to the time reversibility of discrete interactions, as 0j = 
uAtj, but some care must be taken. The generalized discrete trigonometric 
functions are not the truncated series of the standard continuous trigonomet- 
ric functions. They differ in two aspects. Besides being finite, these series 
are defined by the products of, in general, distinct elementary transforma- 
tions. But even when they are equal, as in the particular case discussed in 
the following section, they are not the partial sums of trigonometric series 
which, by the way, are their first order asymptotic limits, as we will see. 
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The Eqs. (10,11,12) imply on 



rt— 1 

cosg 2 9(n) + sing 2 9(n) = ]J (1 + 6)). (22) 

3=0 



CiiGie^r^iiG-^^u^+^-H^Tl ~ smg ^ j ) ( 23 ) 

j=o j=b j=b \sing0(ra) cosg9(n) ) 

and 

712— 1 n2— 1 ni— 1 

n = n n g,_1 (^) = 

i=«i ji=o j2=o 

7 tt i1 / 1 + g)2x-i / cosg6>(^ 2 ) sing#(n 2 ) \ /cosg0(rai) -sing#(m) \ ^ 
j- \— sing#(n 2 ) cosg#(n 2 ) / \sing^(n 1 ) cosg#(rii) /' 

This parallelism with trigonometric functions hides, however, a far more 
larger complexity. According to the definitions Eqs. (17,18), cosg#(n) and 
sing6*(n) are defined by the set {9q, 0i, ■ ■ ■ , n -i}- So, in principle, there is a 
function 9(n), still to be properly defined, but it is already clear that 

9(n 1 )±9(n 2 )^9(n 1 ±n 2 ), (25) 

as 9{n\) and 9(n 2 ) are defined by two distinct sets, both starting with 9 
but ending, respectively with 9 ni -\ and 9 n2 -i. In order to further explore 
the parallelism with the standard continuous trigonometric functions let us 

define, 

cosg[0(rii) =1= 9{n 2 )} := cosg6 l (ni) cosg#(n 2 ) ± sing0(rii) sing#(n 2 ), (26) 
and then, making use of Eqs. (17) and (18), get 

i B2-1 ri2 — 1 m— 1 

cos g [6i(n 1 )-fl(n 2 -i)] = -[n(i+^o+ na-^Kna+^j. ( 27 ) 

j=m j=m ji=o 

i m— 1 ri2 — 1 m— 1 712— 1 

cosg[9(n 1 )+9(n 2 )] = -[ JJ (l+*%) 2 II (1+*%)+ II l 1 "^) 2 II (!-#*)], 

Ji=0 J2=ni ji=0 J2="i 

which reduce, respectively, to Eq. (22) and to 

cosg2#(n) = cosg 2 #(n) — sing 2 #(n) 
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for ni — n\ — n. Then we can see from Eq. (27) and Table 1 that, as a 
particular case of Eq. (25), 

9{n) - 9{n) + 6(0) 

and that 

n-1 

cosg[0(ra) - 0(ra)] = WH + e]). 

3=0 

These peculiarities make the generalized discrete trigonometric functions are 
mathematically interesting consequences of the quasi-unimodularity of the 
evolution function G{6j), Eq. (10) and of the definition (26); they disappear 
if a normalized G{6j) is adopted, but this implies on new equations of motion 
and will be discussed elsewhere. Anyway, these deviations from the standard 
behavior are negligible (cannot be physically detected) in the effectivelly 
continuum-interaction interval (1). 

3 Combinatorial trigonometric functions 

The complete specification of the generalized discrete trigonometric functions 
above depends on the specification of the DHO model. This is realized when 
each time increment Atj is given as a known function of rj. The simplest 
case is the one where Atj is just a constant, it does not depend on j. In 
the other more generic cases the rj-dependence in Atj may be absorbed in 
uj by redefining it into Uj, so that they describe harmonic oscillator with a 
variable, r-dependent, frequency. Let us consider the simplest case and make 

Atj = a, (28) 

where a is a small constant, 

6j = 6 = ua « 1, G(6j) = G{0). 

The function 6{n) just reduces to n6, and the real functions in Eq. (19), the 
generalized discrete trigonometric functions, reduce to the, let us name it, 
discrete combinatorial trigonometric functions 

cose n9±i sine nO = (l±i9) n = ^(±^) J '( n ), (29) 
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(31) 



1 [n/2] (n\ 

coscn^-[(i + ^r + (i-^n = E(-ir^( 2s j ; (30) 

1 [n/2] ( n \ 

sine nO ee -[(1 + ie) n - (1 - = E ("l)^ 2s+1 [ 2fi + J • 

So, according to Eq. (22), they satisfy 

cosc 2 n# + sinc 2 n# = (1 + 6 2 ) n . (32) 

A similarly generalized discrete Bessel function as a finite combinatorial series 
is defined in [5]. Combinatorial trigonometric functions are simpler than 
the generalized discrete trigonometric functions and, therefore, closer to the 
trigonometric functions. As now 9(rii)+9(n 2 ) = 9(rii+n 2 ), it is not necessary 
to define cosg[6'(n 1 ) =F#(n 2 )] and sing[6'(n 1 ) =F#(n 2 )]. The following identities 
are satisfied. 

cose n\6 cose n 2 9 — sine n\ 9 sine n 2 = cose (rii + n 2 )9. (33) 

sine n\9 cose n 2 9 + cose n\9 sine n 2 9 = sine (n x + n 2 )9. (34) 

cose n\9 cose n 2 9 + sine n\9 sine n 2 9 = (1 + # 2 ) ni cosc(n 2 — n 1 )9. (35) 

sine n\9 cose n 2 9 — cose n\9 sine n 2 # = — (1 + # 2 ) ni sinc (n 2 — ni)0, (36) 
and they lead, for n 2 = and rii = n, to 

cosc(-n0) = (1 + fl 2 )- n coscn#, (37) 

sinc(-n0) = -(1 + 9 2 y n sine n9, (38) 
which agrees with Eq. (11). Also 

sinc2n# 

tancn# = — — - — — . (39) 

cosc2n9 + (1 + 9 2 ) n v ; 

These expressions become identical to the ones for the trigonometric func- 
tions if the ubiquotus term (1 + 9 2 ) n is replaced by 1, an anticipative sign of 
the upper limit in Eq. (1) for retrieving the continuum functions, 

(1 + 9 2 ) n « 1 +n9 2 « 1. 
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Table 2: The first five discrete combinatorial trigonometric functions for the 
values of n 



cose 0$ = 


1 




sine 0$ 


= 




cose 19 = 


1 




sine 19 


= 9 




cose 29 = 


1 


-9 2 


sine 29 


= 29 




cose 39 = 


1 


-39 2 


sine 3$ 


= 39 


-9 3 


cose 49 = 


1 


- 6# 2 + 9 A 


sine 49 


= 49 


-49 3 



Table II lists the first five discrete combinatorial trigonometric functions. 

For a discrete hyperbolic harmonic oscillator, one with repulsive interactions 

Ap n ^ = fcf n _iAt n _i, (40) 

which corresponds to an inverted (repulsive) harmonic potential, the solu- 
tions are the generalized iscrete hyperbolic functions 

i n— 1 n— 1 

cosh g ^ = -(n(i+^)+ri(i-^)) (4i) 

3=0 3=0 
i n— 1 n— 1 

smhg£ n = -(11(1 + 9 3 ) - 11(1 - 9,)), (42) 

j=0 j=0 

that are associated to discrete hyperbolic rotations in phase space 



coshg^(n) sinhg^(n) \ 

-smhg0(n) coshg6>(n) / AJL ^ v ; 



They satisfy 



coshg 2 9 n - sinhg 2 9 n = I](l + 0?)- (45) 

j"=o 
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3.1 Combinatorial identities 

Eqs. (30,31,32) or (41,42,45) imply on some useful combinatorial identities. 
Expanding each term in both sides of Eq. (32), 

[n/2] 

cose 2 



2[n/2\ k / \ / \ 

^^r ,w t-VJt + ) <47) 

{i+0 2 r = p{^0 2 \ (48) 

and equating the coefficients of the 9 2k terms in both sides of Eq. (32), we 
have 

H)'EL"J(:]-iL.v,)f",))=(i]. («> 



t=0 



2(k-t))\2t) ^ \2(fc-t) - 1/ V2t+ 1/ 



where, there is a clear separation, in left hand side, of even and odd terms, 
a trace of sines and cosines. It can be put in a more symmetrical form as 

3.2 Retrieving the continuum limit 

For completeness sake we make here a brief review [1] of how the effective 
continuum limit is reached. For large values of n the following approximation 
is valid 



k! ~ kV f ° r n>>1 ' (51) 



and this turns Eq.(29) into 



— 77 

cose nd±i sine nd = (l±id) n « ^(±i0) J — , (52) 

3=0 
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and then the Eqs.(30) and (31) become, respectively 



[n/2] n 2s 



coscnO ks ^2 (— 1) s («cj) 2s ^ ^ cos(acun), (53) 

K2] n 2s+l 

sincn# ps (-l) s ((acj)) 2s+1 — — - ps sin(acjn). (54) 

The middle terms of Eqs. (53) and (54) are partial sums of the series repre- 
senting the respective trigonometric functions which are then the asymptotic 
limits of the respective largest contributions (Eq. (51)) from each term of 
the combinatoric function series for n — > oo. The continuous variable time 
enters in the description of the motion as a discrete parameter (as it happens 
also to the position variable) 

t n = t + na. (55) 

It is important to underline that the effective standard picture of a continu- 
ously interacting harmonic oscillator emerges, in a natural and unavoidable 
way, as the number n of interactions grows, irrespective of the fixed values of 
lu and of a — At. The smallness of At is never required. That is what char- 
acterizes this discrete-to-continuum transition as being of the kind described 
by the Bohr Correspondence Principle. Just for the sake of comparison we 
remind well known standard limiting procedures like, for example, to replace 
in Eq. (29) At by ^ which, naturally goes to zero as N grows unboundly: 

cose coAt±i sine uj At = (l±i^)"xe iuAl . (56) 

N, in contradistinction to n, is not the number of interactions, which is 
meaningless for continuum interaction; it is just an ancillary mathematical 
artifact. Such limiting procedures with time lapse approaching zero with the 
growth of N, lead to a genuine continuum interaction and the last term in 
the right hand side of Eq. (56) represents an exact, not just an asymptotic, 
expression. For discrete interactions, in contradistinction, where At is fixed, 
it is just a first order approximation, and there will always be known lower 
order residues. The larger is n the better the first order contribution (in the 
sense of Eq. (51)) series approach their respective trigonometric functions. 
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As 6 << 1, the second order contributions are always negligible face the 
first order ones, but for large enough values of n they may become physically 
detectable producing new physical effects not contained in the continuum 
formalism. 

4 Discrete angular momentum and Hamilto- 
nian 

— * 

For the continuum harmonic oscillator the angular momentum L := r x p 
is time independent, and the generator of infinitesimal rotations. For the 
discrete harmonic oscillator, with the quasi- unimodular elementary trans- 
formations (10), the above standard definition for the angular momentum 
makes it explicitly time dependent (6 = uAt) 

n-1 

L n ■= f n xp n = (l + = II (! + 9 !)^o, (57) 

3=0 

and, besides, since there is no infinitesimal rotations (only polygonal trajec- 
tories) it cannot be their generators. Similarly, the Hamiltonian, i.e. the 
kinetic energy plus the potential energy, is also a constant of motion and 
a generator of infinitesimal translation in time for the continuum harmonic 
oscillator. For the discrete one, however, there is no potential energy, only 
kinetic energy, and its time evolution is trivial (free evolution) except on a 
finite set of interaction points. So, the meaning of the concept of a Hamil- 
tonian in a discrete interaction context must be reassessed. Let us consider 
the discrete Hamiltonian 

\ 2 -r 2 J, 

written in terms of the homogenized units of Eq. (7). According to Eqs. (8), 
(9) and (22), 

H n =J[(l + e])H Q , (58) 

3=0 

and so, H n is also explicitly time-dependent (9j = wAtj). Actually, H n , 
and also L n , would be more properly described as being n, rather than t n , 
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dependent. They remain constant in the time lapse At n between the n — th 
and the (n + 1) — th interactions. At each interaction they get larger by a 
factor 1 + 0*. 

The discrete time evolution (10) does not preserve the canonical relations 
either. The Poisson brackets are not preserved 

where F n and G n are functions of r n and p n . For example, 

n-1 

{Pn,r n } = {po,r }U( 1 + d ")- (6°) 

j=0 

— * 

This non preservation of H n , L n and of the Poisson bracket is a consequence 
of the non-unimodular character (10) of the elementary transformation. 

Whereas the condition n » 1, as an inferior bound, is required for 
retrieving the standard continuous trigonometric functions as first order ap- 
proximations for the above DHO solutions, the preservation of the canonical 

— # 

relations, or equivalently, the conservation of H n (and of L n ) imposes the 
upper bound condition 

nB 2 « 1, 

of the relation (1) for the validity of continuum interaction as a good first 
order approximative description. Eq. (57), written as 

— * — * 

- ~ ^ = ^At n ^Ln-l (61) 

shows that if Ai n _! could be made infinitely close to zero we could have 
a time derivative (if the ratio on the left hand side remains well defined) 
and L n would be conserved. The same goes for H n . But for the discrete 
oscillator Atj is fixed, it cannot go to zero. The emerging potential energy, 
as an effective concept, corresponds to the energy of the still travelling (not 
yet absorbed) exchanged quanta. 

5 Physical interpretations and Conclusions 

It is plausible that in our quotidian Atj, the time interval between two con- 
secutive discrete interactions, being extremely small, uiAt n « 1, be beyond 
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any todays conceivable experimental detection. Therefore, even if there is 
no legitimate infinitesimal transformation but only finite ones, in a practical 
sense of approximation, at the measure of uAt n « 1, they can be treated as 
if they were infinitesimal [1]. Rigorously H n and L n , with quasi-unimodular 
elementary transformations, are not time-independent, but just quasi. 

H n — (1 + 9j)H n -i PS -f/n-l 
— » 2 - * - * 

L n = (1 + ^ n _ 1 )L n _i ~ £ n -i- 

Up to a first order approximation both, H n and L n appear to be time 
independent and the system can be taken as a continuum one. 

We list, after all these considerations, three possible alternative scenarios: 

1. The time evolution equation (8) is correct but Eqs. (57) and (58) are 
not the correct expressions for, respectively the angular momentum 
and the energy of a DHO. They should have a factor 1/ ydTj=o (1 + @f) 

that would make both H n and L n , time independent but still could not 
preserve the canonical relations. The preservation of the canonical re- 
lations with all its implications (vide quantum mechanics, for instance) 
would be just a first order approximation artifact and not a funda- 
mental truth of nature. This would be, of course, a very controversial 
interpretation; 

2. The equations (8), (57) and (58) are correct. The explicit time depen- 
dence of Eqs. (57) and (58) shows that they do not describe a closed 
system. An energy and angular momentum supplier is missing in this 
description. On the other hand, all the terms in these equations are 
definite positive which means that is it not possible to build a closed 
(energy and angular momentum conserving) system only with discrete 
classical harmonic interactions. The concept of energy-and-angular- 
momentum conserving classical harmonic oscillator, rigorously, does 
not exist in a discrete interaction context; it is just another asymptotic 
physical idealization like many others (ideal gas, reversible process, 
etc.) in physics. 

3. Equation (8) is not correct; it needs to be normalized into 

VPn/ Jl + Ql_ x V -0 n -l 1 J \Pn-lJ 
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Then the new elementary transformation 

G(gn-l)= / 1 ( / V) (63) 

is unimodular, leads to normalized generalized discrete trigonometric 
functions, with 

cosg 2 6 {n) + sing 2 9{n) — 1, (64) 

replacing (22), and so, preserves the Hamiltoninan, the angular mo- 
mentum and the Poisson Bracket structure. We do, certainly, favor 
this alternative but it leads to new second order (on 9) effects that, 
in principle, can and should be experimentally checked. Its discussion 
will be presented elsewhere in a quantum mechanics context. 

A fourth alternative, of course, could be quantum exchange being just an 
approximation for an actually continuously interacting world. As it is contra- 
dicted by all the great theoretical and experimental achievements of quantum 
theory it was not even listed. In summary, the possible idea of a fundamen- 
tally discrete interaction puts in check the entire structure of physics based on 
interaction continuity — is it exact or just an approximation experimentally 
validated by the smallness of the interaction quanta? 
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